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NUMERICAL  SIMULATION  OF  CONFINED 
MULTIPLE  TRANSVERSE  JETS 


Farhad  Davoudzadeh^ 
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ABSTRACT 

Behavior  of  unconfined  transverse  jets  has  been  studied  extensively,  but  little  work  is 
reported  on  the  flow  characteristics  of  confined  transverse  jets.  The  latter  has  been  numerically 
investigated  using  a  number  of  RANS  codes.  The  computational  results  obtained  from  these  codes 
have  been  evaluated  against  the  existing  experimental  data,  and  an  LES  results  reported  in  the 
literature.  Furthermore,  an  extensive  validation  effort  has  been  conducted  to  characterize  the 
performance  the  codes  for  predicting  the  flow  within  a  propulsion-related  mixing  configuration.  The 
validation  case  involves  eight  circumferentially  spaced  transverse  jets  issuing  radially  into  an 
axisymmetric  main  flow,  a  configuration  relevant  for  gas  turbine  burners  and  new  liquid  rocket 
engine  preburners.  The  main  flow  Reynolds  number  was  1.7  x  10^  and  the  jet-to-main  flow 
momentum  flux  ratio  was  sixteen.  The  momentum  and  scalar  mixing  was  investigated  through  the 
solution  of  the  Reynolds- Averaged  Navier  Stokes  (RANS)  equations.  The  solutions  of  three 
commercial  RANS  solvers.  Fluent,  STAR-CCM+,  and  CFD++,  are  compared  to  experimental  data 
and  large-eddy  simulation  (LES)  results  available  in  literature.  Due  to  demonstrated  periodicity, 
only  a  one-eighth  pie-shaped  section  of  the  geometry  was  considered.  The  different  commercial 
codes  used  the  same  geometry,  grid,  boundary  conditions,  and  variations  of  the  k-e  turbulence 
model.  The  LES  results  obtained  from  literature  used  a  different  grid,  but  the  same  geometry.  All 
numerical  simulations  using  the  above  mentioned  codes  capture  salient  flow  structures  such  as  the 
counter-rotating  vortex  pair  (CRVP).  Experimental  data  used  for  validation  of  the  codes  include 
mean  axial  velocity  and  jet  fluid  mixture  fraction  profiles  (at  three  distinct  axial  locations),  jet 
trajectory,  turbulent  kinetic  energy  distributions,  and  velocity  and  mixture  fraction  cross-plane 
distributions.  All  CFD  results  except  CFD++,  exhibit  symmetrical  solutions  about  the  center  plane. 
The  current  investigation  shows  that  although  all  codes  considered  predict  the  experimental  data 
with  various  degrees  of  accuracy.  Fluent  using  the  standard  k-e  turbulence  model  with  the  standard 
wall  function,  and  LES  results  compare  exceptionally  well  with  the  experimental  data  for  this  flow 
regime  and  configuration. 


Nomenclature 


A,  b,  c 

empirical  constants 

C 

local  mean  mass  fraction  of  the  jet 

C 

average  mass  fraction  of  the  jet  over  the  cross  sectional  area  A 

D 

diameter 

J 

momentum  flux  ratio 

m 

mass  flow  rate 

n 

number  of  jets 

Q 

volume  flow  rate 

r 

velocity  ratio 

X 

streamwise  coordinate 

y 

cross  stream  coordinate 

j 

subscript  for  jet  flow 

Us 

Unmixedness 
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I.  Introduction 


Transverse  jets,  or  jets  in  cross  flow  (JICF),  provide  an  effective  mixing  mechanism,  via  vortical 
structures,  that  can  be  utilized  in  a  number  of  applications.  Examples  of  these  applications  include  mixing 
of  diluents  with  core  fluid  in  a  liquid  rocket  engine  prebumer  or  gas  turbine,  vertical  and  short  take-off  and 
landing  (V/STOL)  aircraft,  and  film  cooling  of  turbine  blades.  Past  studies  have  demonstrated  that 
excessive  temperature  variations  at  the  turbine  inlet  have  a  detrimental  effect  on  turbine  life  and 
performance.  In  staged  combustion  rocket  engines,  the  function  of  a  prebumer  is  to  provide  combustion 
gases  that  will  drive  the  turbine,  which  in  turn  drive  the  pumps,  before  being  exhausted  into  the  main 
combustion  chamber.  Survival  and  durability  of  the  turbine  blades  depends  upon  thermal  stresses  applied 
by  the  combustion  hot  gases  exhausted  by  the  preburner.  The  cooler  and  the  more  uniform  the  combustion 
gases  the  better  the  survivability  of  the  turbine  blades.  In  this  application  the  LOX  diluent  jet  in  the  cross 
flow  can  provide  enhanced  mixing  mechanism  of  the  combustion  gases  to  provide  a  uniform  inlet  to  the 
turbine. 


The  rocket  preburner  offers  a  number  of  technical  challenges  for  CFD,  including  turbulent 
combustion,  complex  geometry  with  strong  streamline  curvature,  very  large  density  ratio  between  the 
mixing  fluids,  and  transcritical  phenomena  due  to  the  high  operating  pressure.  The  current  validation  effort 
is  essentially  the  first  phase  of  an  overarching  process  of  validating  codes  for  rocket  preburners,  and 
employs  currently  available  experimental  data  associated  with  a  gas  turbine  burner  configuration  [1,  2]. The 
experimental  setup  involved  an  eight  JICF  configuration  issuing  into  a  fully-developed  pipe  flow  with  both 
fluids  being  air  nominally  at  standard  temperature  and  pressure.  The  purpose  of  this  phase  of  the  validation 
effort  is  to  establish  a  foundational  methodology  that  is  able  to  accurately  simulate  the  complex  geometry 
in  a  preburner-like  configuration.  This  will  essentially  establish  a  departure  point  as  we  incrementally  add 
salient  features  of  the  preburner,  including  density  ratio,  high  pressure,  and  chemistry. 

The  experimental  work  was  numerically  modeled  using  three  different  computational  fluid 
dynamics  Codes  (CFD):  Fluent,  CFD++,  and  STAR-CCM+.  The  computational  work  was  performed  to 
identify  the  current  capabilities  of  the  CFD  tools  in  modeling  the  complex  flow  structures  encountered  in 
JICF,  and  to  identify  the  appropriate  CFD  model  elements  (solver,  turbulence  model,  computational  extend 
and  grid,  etc.)  related  to  the  flow  regime  under  consideration.  Current  simulations  show  a  good  comparison 
between  the  computational  results  obtained  by  Fluent  and  STAR-CCM+  and  identifies  relevant  CFD 
parameters  (turbulence  model,  geometry  and  grid  distribution,  solver  type)  of  importance  in  this  flow 
regime.  Efforts  are  continuing  to  perform  more  detailed  calculations  looking  into  grid  resolution,  boundary 
conditions,  and  solver  models. 

Demuren  [3]  categorized  four  different  models,  in  terms  of  accuracy,  that  has  been  used  in 
investigations  of  the  JICF:  empirical,  integral,  perturbation  and  numerical.  The  empirical  method  is  the 
simplest  model  that  relies  on  the  correlation  of  experimental  data,  and  the  accuracy  of  an  empirical  model 
will  in  general  depend  on  whether  a  given  point  of  interest  lies  within  the  cloud  of  data  used  for 
construction  of  the  correlation;  the  empirical  model  is  essentially  a  curve  fit.  Due  to  their  low  cost  and  ease 
of  use,  empirical  models  are  most  useful  for  first-order  estimates  and  as  qualitative  checks  for  results 
produced  by  other  methods. 

One  of  the  the  most  common  JICF  attributes  treated  using  empirical  models  is  the  jet  trajectory. 
For  a  single  circular  turbulent  jet  injected  normally  into  a  cross  flow,  the  trajectory  has  the  form  [3]: 


f  ^ 


D 


T 

j 


where  D  is  the  jet  diameter,  J  is  the  momentum  flux  ratio,  x  and  y  are  the  cross-stream  and  streamwise 
coordinates  of  the  jet  trajectory,  and  a,  b,  and  c  are  empirical  constants.  In  the  range  of  J  between  2  and 
2,000,  a  has  a  value  between  0.7  and  1.3,  b  has  a  valuebetween  0.36  and  0.52  and  c  takes  a  value  between 
0.28  and  0.40,  depending  on  experimental  conditions  [3]. 
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On  the  other  end  of  the  spectrum,  Demuren  [3]  pointed  out  that  “numerical  models  have  the  most 
potential  for  wide  generality  and  can,  in  principle,  be  applied  to  the  whole  range  of  jet  in  cross  flow 
situations,  confined  or  unconfmed,  low  medium  or  high  J,  single  or  multiple  jets,  impinging  on  a  wall  or 
on  other  jets,  swirling,  homogeneous  or  heterogeneous  cross  flow,  compressible  or  incompressible,  etc.” 

The  research  results  presented  in  this  paper  numerically  investigates  the  various  flow 
characteristics  experienced  when  jets  and  cross  flow  interact,  studies  the  induced  resultant  flow  structures 
and  their  contribution  to  the  mixing  of  the  propellants.  For  validation  and  verification  of  the  simulations, 
the  numerical  work  uses  several  RANS  codes  and  compares  the  simulation  results  against  experimental 
data  and  a  large-eddy  simulation  (LES)  available  in  the  literature. 

II  Geometry  and  Grid 

The  experimental  data  used  for  validation  of  the  computational  work  were  obtained  from  the  work 
performed  by  A.  Strzelecki,  et  al.  [2-3].  In  these  experiments,  the  mixing  of  eight  isothermal  jets  issuing 
into  a  circular  pipe  flow  is  investigated  by  means  of  laser  Doppler  anemometry,  particle  image  velocimetry, 
and  planar  laser-induced  fluorescence  techniques. The  measured  variables  include  velocity  and  scalar 
concentration  fields.  Computational  representation  of  the  geometry,  in  the  vicinity  of  the  diluent  injection 
ports,  is  shown  in  Figure  1. 


Diluents  Flow 


Figure  1.  Geometry  in  the  vicinity  of  the  jet  injection  system. 

It  should  be  noted  that  this  gemetry  is  constructed  for  the  CFD  analysis  and  is  considered  to 
closely  represent  the  experimental  configuration.  A  similar  simplification  of  the  geometry  is  also  adapted 
by  the  original  investigators  [3]  for  their  LES  calculations. 

Figure  2,  shows  the  grid  distribution  (a)  on  and  around  an  injection  port,  and  (b)  at  the  inlet/exit  of 
the  core  flow.  A  total  of  approximately  800k  hexahedral  only  elements  are  used.  As  shown  in  the  inlet/exit 
grid  distribution,  singularity  axis  is  used  in  the  construction  of  the  grid.  Wall  grid  distribution  corresponds 
to  ay+  of  0.71  -1.2. 
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Figure  2.  (a)  Grid  distribution  on  an  injection  port  (b)  and  at  the  core-flow  inlet/exit 


COMPUTATIONAL  DOMAIN  AND  BONUDARY  CONDITIONS 

Figure  3  shows  the  computational  domain  and  the  boundary  conditions  used  in  the  simulation.  To 
take  advantage  of  symmetry  in  the  geometry,  only  a  forty-five  degree  pie  section  of  the  geometry  is 
simulated,  with  the  planes  of  symmetry  specified  at  each  side  of  the  pie,  as  shown  in  the  figure.  The 
upstream  boundary  is  located  lOD  upstream  of  the  diluent  injection  port,  where  D  is  the  main  pipe 
diameter.  The  downstream  boundary  condition  is  similarly  located  lOD  downstream  of  the  injection  port. 
The  upstream  BC  is  specified  at  an  inlet  velocity  of  25  m/s.  The  inlet  jet  velocity  is  58.1  m/s.  Considering 
the  contraction  of  the  jet  tube  before  intersection  with  the  main  pipe,  the  jet  velocity  at  the  jet  orifice  is  100 
m/s.  This  dictates  a  velocity  ratio(jet  to  main  flow)  to  be  (100/25=)  4,  that  corresponds  to  a  momentum  flux 
ratio  of  16.  The  exit  BC  uses  the  static  pressure  and  is  set  to  standard  atmospheric  pressure.  The  main  core 
diameter  is  100  mm,  the  jet  diameter  is  8.1  mm  that  contracts  to  6.1  mm  at  the  orifice.  These  conditions 
match  that  for  the  experimental  and  LES  data.  The  fluid  used  for  both  diluent  inlet  and  for  the  main  core 
flow  is  air 


Diluents'  Inlet 


Figure  3.  Computational  domain  and  boundary  conditions. 


Ill  Flow  Solver  and  Operating  Conditions 

All  simulations  are  3D,  second-order  accurate,  using  Ansys-Fluent  (pressure-velocity 
coupled)  flow  solver  Release  13.0.,  CFD++,  AND  STAR-CCM+.  In  addition,  the  LES  results  from  Priere 
et  al.  [3]  are  also  considered.  Equations  solved  are  three  for  momentum,  energy,  turbulent  kinetic  energy 
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(TKE),  turbulent  dissipation  rate  (TDR),  species,  and  continuity.  The  second-order  spatial  discretization 
solution  method  was  used  in  all  calculations. 

The  inlet  boundary  conditions  (BC)  specify  velocity,  temperature,  and  turbulence  intensity.  The 
exit  BC  sets  static  pressure  and  temperature.  Walls  are  no  slip  and  adiabatic.  Turbulence  models  considered 
include  the  standard  and  realizable  two  equation  k-e  models  using  standard,  enhanced,  or  the  low  Reynolds 
number  wall  treatment  options. 

All  the  results  presented  are  after  the  numerical  iteration  has  converged.  Convergence  is 
considered  achieved  when  the  residuals  of  the  equations  being  solved  have  dropped  by  approximately  four 
orders  of  magnitude,  gradient  of  the  residual  curves  had  approached  zero  and  have  plateaued,  there  are 
negligible  changes  in  the  flow  variables,  and  mass  is  conserved  throughout  the  computational  domain.  All 
calculations  are  performed  on  a  four  processor  workstation.  Simulations  typically  took  approximately 
10000  iterations  to  achieve  a  fully-converged  solution.  This  typically  took  about  24  hours,  dependent  on 
how  many  grid  points  were  used.  On  quad-core  Intel  Zeon  processor  the  clock  time  used  was  25  x  10-3 
microseconds  per  iteration  per  grid  point. 


IV  Computational  Results  and  Experimental  Data 
FLOW  CHARACTERISTICS: 

Flow  characteristics  of  the  JICF  are  described  in  detail  by  Fric  &  Roshko  [4].  The  flow  is 
dominated  essentially  by  four  different  types  of  vortical  structures  as  shown  in  figure  4:  (1)  a  counter 
rotating  vortex  pair  (CRVP).  This  vortex  pair  is  the  most  dominant  flow  feature  in  the  JICF.  It  is  generated 
by  the  realignment  of  jet  shear  layer  vorticity  through  a  complex  interaction  with  the  cross  flow  [5].  The 
CRVP  is  the  primary  flow  structure  responsible  for  mixing  of  the  jet  and  the  core  flow  (2)  Horseshoe 
vortex:  generated  at  the  root  of  the  jet  as  it  enters  the  cross  flow.  This  type  of  vortex  structures  is  a  typical 
flow  pattern  that  is  created  around  a  bluff  body  mounted  normal  to  a  surface.  It  is  generated  due  to  the 
adverse  pressure  gradient  on  the  windward  side  of  the  jet  and  the  adjacent  wall  due  to  the  dynamic  pressure 
profile  within  the  boundary  layer.  (3)  Jet  shear  layer  vortices:  these  vortices  originate  from  instability  of  the 
separated  boundary  layer  at  the  orifice  of  the  jet  and  evolve  around  the  surface  of  the  jet.  The  shearing  of 
the  core  flow  via  the  jet  flow  generates  Kelvin-Helmholtz  instability  causing  the  flow,  in  the  shear  region, 
to  roll  up  to  a  ring  vortex  around  the  jet.  (4)  Wake  vortices:  these  vortices  are  found  in  the  wake  of  the  jet 
and  they  connect  the  CRVP  to  the  wall.  These  structures  where  initially  investigated  experimentally  by  Fric 
and  Roshko  [4]. 


Counter  Rotaiing 
Vortex  Pair 


Figure  4.  Schematic  of  various  vortex  structures  in  JICF,  After  Fric  &  Roshko  [3] 

The  RANS  numerical  simulation  using  ANSYS-Fluent  captured  some  of  these  vortical  structures 
as  described  in  the  following.  Figure  5  shows  an  isometric  view  of  the  flow  field.  The  characteristic 
velocity  Uo  is  the  mean  velocity  at  the  main  pipe  inlet.  The  non-dimensional  axial  velocity  gives  an 
indication  of  the  extent  of  the  jet.  The  figure  also  shows  three-dimensional  streamlines  in  the  cross- 
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sectional  planes  normal  to  the  axial  direction.  These  streamlines  delineate  the  CRVP.  Expansion  of  the 
CRVP  in  the  axial  direction  is  demonstrated  in  these  cross-sectional  streamlines. 

Diluent  Flow 


♦ 


'  e- 


Figure  5.  Axial  velocity  &  cross  sections  including  streamlines  showing  CRVPs’  structure 

(ANSYS/Fluent  simulation) 

VELOCITY  DISTRIBUTION: 

Figure  6  shows  contours  of  streamwise  velocity  at  axial  cross  sections  2d,  5d,  and  lOd  (d  is  the 
diameter  of  the  jet)  downstream  of  the  jet  for  the  computational  results  obtained  using  Fluent,  STAR- 
CCM+,  CFD++,  FES,  and  the  experimental  data.  Comparing  the  velocity  distributions,  it  is  evident  that 
FES  and  Fluent  results  compare  very  well  with  the  experimental  data  for  all  the  axial  locations  considered. 
The  RANS  simulations  shown  in  the  figure  were  run  using  the  standard  k-e  model  turbulence  model  with 
standard  wall  function,  except  CFD++  where  the  standard  k-e  model  was  not  available  and  the  realizable  k- 
G  model  with  the  wall  function  was  used.  A  noticeable  feature  in  the  CFD++  results  is  the  asymmetry  in  the 
solution  that  becomes  very  pronounced  at  the  axial  station  of  x  =  1  Od.  This  was  brought  to  the  attention  of 
Metacomp  Technologies,  the  developers  of  the  CFD++  software,  in  a  private  communication. 
Consequently,  Metacomp  ran  a  calculation  in  which  they  indicated  that  “in  order  to  obtain  a  non-oscillating 
steady-state  solution,  flux  blending  was  used  with  75%  second-order  flux  blended  with  25%  first-order 
flux”.  Metacomp  also  indicated  that  when  run  in  higher-order  transient  mode,  this  type  of  flow 
configuration  has  a  tendency  to  develop  asymmetries  which  should  be  expected  to  change  the  time- 
averaged  interior  flow- field,  including  the  exit  mass  fluxes. 
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Figure  6.  Axial  velocity  distribution  at  different  axial  stations 

Figure  7  compares  computational  radial  profiles  of  axial  velocity,  obtained  via  the  above 
mentioned  codes,  with  the  experimental  data  at  axial  locations  of  x  =  2d,  5d,  and  lOd  downstream  of  the 
injection  port,  where  d  is  the  jet  orifice  diameter.  At  the  axial  station  of  x  =  2d,  all  codes  predict  the 
experimental  profiles  of  the  axial  velocity  very  well,  but  for  the  locations  of  x  =  5d,  and  x  =  lOd  the  LES 
and  the  Fluent  predictions  are  much  closer  to  the  experimental  data  than  CFD++  and  STAR-CCM+. 
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Figure  7.  Predicted  vs.  experimental  radial  profiles  of  the 
axial  velocity  at  different  axial  locations 


Figure  8  shows  comparison  of  the  experimental  data  with  the  Fluent  predictions  using  the 
standard  k-e  and  the  realizable  k-e  models,  and  the  CFD++  predictions  using  the  realizable  k-e  model.  The 
figure  demonstrates  that  the  Fluent  and  the  CFD++  predictions  show  similar  patterns  when  the  realizable  k- 
G  model  is  used.  This  could  be  attributable  to  the  turbulence  models  constants  used  in  the  standard  k-e  and 
the  realizable  k-e  models.  The  standard  k-e  models  uses  five  constants,  namely  Cp  =  0.09,  Cle  =  1.44,  C2g 
=  1.92,  ak  =  1.0,  and  ae  =  1.3  whereas  in  the  realizable  k-e  model  Cle  and  Cp  are  not  constants  [6].  In 
particular,  the  maximum  value  that  can  be  allowed  in  the  realizable  k-e  model  for  Cp,  in  the  whole 
computational  domain,  is  0.09.  In  most  cells  this  value  is  smaller  than  0.09.  Considering  the  transport 
equations  for  the  k-e  turbulence  model,  the  higher  value  of  the  Cp  would  lead  to  a  higher  production  of  the 
turbulent  kinetic  energy  when  the  standard  k-e  is  used  versus  the  realizable  k-e  (see  Figure  9),  and  this  may 
be  responsible  for  the  better  predictions  when  the  standard  k-e  model  is  used. 
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Figure  8.  Experimental  data  vs.  predicted  radial  profiles  of  the  axial  velocity  using  standard  k-e  and 

realizable  k-e 


(a)  (b) 


Figure  9.  Predicted  turbulent  kinetic  energy  with  Fluent  (a)  Standard  k-e  with  standard  wall 
function,  (b)  realizable  k-e  with  enhanced  wall  treatment 
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There  is  a  disagreement  between  a  significant  number  of  authors  about  two-equation  turbulence 
models  and  their  application  for  different  flow  regimes. 

As  Menter  [7]  has  pointed  out,  “there  is  no  agreement  on  the  standards  by  which  to  measure  the 
improvements  achieved  by  proposed  new  models,  or  alterations  to  the  existing  models.  Many  times  new 
models  are  based  on  theoretical  concepts,  which  by  themselves  involve  severe  assumptions  about  the  nature 
of  turbulence.  It  is  often  unclear  whether  the  improvements  presented  for  one  type  of  flow  (e.g.,  boundary- 
layer  flows)  will  not  lead  to  deterioration  for  another  class  of  equally  important  flows  (e.g.,  free  shear 
flows).”  For  example,  an  argument  had  been  made  by  Pope  [8]  in  1978  that  the  standard  k-e  model 
overestimates  spreading  of  the  round  jet  and  suggested  corrections  to  the  model.  However,  as  pointed  out 
by  Rubel  [9],  the  Pope  correction  has  an  adverse  effect  on  the  model  predictions  for  the  radial  jet  and 
underestimates  the  radial  jet  spreading  by  60%. 

Again  as  Menter  [7]  has  stated  “none  of  the  available  theoretical  tools  (dimensional  analysis, 
asymptotic  expansion  theory,  use  of  direct  numerical  simulations  (DNS)  data,  renormalization  group 
(RNG)  theory,  rapid  distortion  theory,  etc.)  can  provide  constants  to  that  degree  of  accuracy. 

The  only  way  to  establish  the  validity  of  theoretical  arguments  is  to  carefully  test  the  resulting 
model  against  a  number  of  challenging  and  well-documented  research  flows.” 

The  comparison  with  the  experimental  data  is  usually  believed  to  allow  some  conclusions  about  a 
model's  ability  to  perform  in  engineering  applications. 


SCALAR  CONCENTRATIONS: 

Figure  10  shows  the  jet  species  mass  fraction  distribution  at  the  longitudinal  cross  section  for  the 
experimental  data  and  the  computational  codes  considered.  The  experimental  data  shows  a  more  rapid 
expansion  of  the  jet  scalar  field  than  STAR-CCM+  and  CFD++.  Fluent  results  again  show  a  better 
matching  with  the  experimental  data  than  the  other  two  codes. 


Experiment  Fluent  Star-ccm+  CFD'^'^ 


Figure  8.  Jet  mass  fraction  distribution 


Experimental  and  predicted  concentration  radial  profiles  at  the  axial  locations  of  x  =  2d, 
5d,  and  lOd  are  shown  in  Figure  11  that  similarly  shows  the  Fluent  results  are  in  much  better  agreement 
with  the  experimental  data  than  STAR-CCM+  or  CFD++. 
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Figure  9.  Computational  and  experimental  jet  concentration 
profiles  at  x  =  2d,  5d,  and  lOd 


Figure  11  shows  the  degree  of  mixing  of  the  jet  flow  with  the  core  flow.  The  unmixedness 
parameter  Us  is  defined  as: 

2  -[(C-CjdA 

"  c(i-c)  c(i-c) 

where  ac  is  the  statistical  variance  of  the  jet’s  mass  fraction  in  the  cross  plane  at  a  given  axial 

location  (X/d),  C  is  the  local  mean  mass  fraction  of  the  jet,  and  C  is  the  average  mass  fraction  of  the  jet 
over  the  cross  sectional  area  A. 

Consistent  with  the  contour  plots  of  the  jet  mass  fraction  shown  in  Figure  10,  Fluent  results 
show  a  higher  degree  of  mixedness  than  the  other  codes. 
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Figure  10.  Experimental  and  predicted  unmixedness 


V  Conclusions 

A  series  of  numerical  simulations  is  performed  to  predict  the  flow  induced  by  eight  isothermal  jets 
issuing  in  a  circular  pipe  flow.  The  computational  codes  used  for  these  simulations  included  an  LES  code, 
Fluent,  STAR-CCM+,  and  CFD++  codes.  The  computational  parameters  were  chosen  to  be  identical  for 
the  codes  under  consideration,  noting  that  some  limitations  were  present  in  the  options  for  the  different 
codes.  The  turbulence  model  used  in  these  simulations  was  the  standard  k-e  model  with  the  standard  wall 
function  for  Fluent  and  STAR-CCM+,  and  the  realizable  k-e  model  with  the  standard  wall  function  for  the 
CFD++;  standard  k-e  is  not  available  in  the  CFD++  code.  Further,  the  computational  results  were  compared 
with  the  experimental  data  for  the  velocity  distributions,  concentration  distributions,  velocity  profiles,  and 
scalar  concentration  radial  profiles  at  various  axial  stations. 

All  simulations  captured  the  dominant  counter  rotating  vortex  pair  (CRVP)  flow  structure,  known 
to  be  present  in  the  JICF  flow  regime.  Quantitative  comparison  of  the  experimental  data  with  the 
computational  results  for  the  radial  profiles  of  the  axial  velocity  and  the  jet  concentration  showed  that 
Fluent  using  the  standard  k-e  predictions  were  closer  to  the  experimental  data  than  either  STAR-CCM+ 
orCFD++.  This  conclusion  is  consistent  with  similar  conclusions  reached  by  Karvinen  et  al.  [10],  involving 
thesimulation  of  a  circular  jet  in  a  cross  flow  using  Fuent  employing  various  turbulence  models.  Karvinen 
concluded  that  results  from  standard  k-e  model  were  closest  to  the  measured  values.  The  turbulence  models 
that  were  evaluated  included  three  high-Reynolds  number  k-e  models,  two  k-co  models,  six  low-Reynolds 
number  k-e  models  as  well  as  the  Reynolds  stress  model.  Menzies  [11]  also  evaluated  various  turbulence 
models  for  a  similar  flow  regime  for  the  isothermal  flow  in  a  gas  turbine  combustion  system  and  concluded 
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that  of  the  turbulence  models  used,  the  realizable  model  performed  poorly,  while  the  RNG  model  did  not 
deliver  any  benefits  to  the  velocity  predictions,  but  the  standard  K-e  model  delivered  acceptable  results  in 
comparison  with  the  measurements. 
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